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Abstract 

We study the effects of finite proton mass on the energy levels of hydrogen 
atoms moving transverse to a superstrong magnetic field B with generalized 
pseudomomentum K±. Field strengths of order B ^ 10^^ Gauss are typically 
found on the surfaces of neutron stars, but we also study the regime B ^ 
Bcrit = 4.23 X 10^^ Gauss, where the Landau excitation energy of the proton 
is large. We adopt two different approaches to the two-body problem in 
strong magnetic field, and obtain an approximate, but complete solution of the 
atomic energy as a function of B and K±. We show that, for B » Bcrit ^ there 
is an orthogonal set of bound states, which do not have any Landau excitation 
contribution in their energies. The states with very large K± have small 
binding energy and small transverse velocity, but are nevertheless distinct 
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from the fully ionized states. The final results for the excitation energies are 
given in the form of analytical fitting formulae. 

The generalized Saha equation for the ionization-recombination equilib- 
rium of hydrogen gas in the presence of a superstrong magnetic field is then 
derived. Although the maximum transverse velocity of a bound atom de- 
creases as B increases, the statistical weight due to transverse motion is ac- 
tually increased by the strong magnetic field. For astrophysically interesting 
case of relatively low density and temperature, we obtain analytic approxi- 
mations for the partition functions. The highly excited bound states have a 
smaller statistical weight than the fully ionized component. 

32.60.-Fi, 97.10.Ld, 31.20.Di, 97.60.Jd 
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I. INTRODUCTION 



It is well known that a superstrong magnetic field of order B ^ 10^^ Gauss, typically 
found on the surfaces of neutron stars, can dramatically change the structure of neutral 
atoms and other bound states (see, e.g., Ref. for an early review and Ref. [0 for a recent 
text on atoms in strong magnetic fields). The atomic unit for the magnetic field strength 
and a dimensionless parameter b are 

^ ^ 2.35 X 10^ G; b=—. (1.1) 

h Bo 

When b >> 1, the cyclotron energy of the electron hue = h{eB/mec) = 11.58-Bi2 keV, where 
Bi2 is the magnetic field strength in units of 10^^ G, is much larger than the typical Coulomb 
energy, thus the Coulomb forces act as a perturbation to the magnetic forces on the electrons, 
and at most temperatures the electrons settle into the ground Landau level. Because of the 
extreme confinement of electrons in the transverse direction, the Coulomb force becomes 
much more effective for binding electrons in the parallel direction. The atom has a cigar-like 
structure. Moreover, it is possible for these elongated atoms to form molecular chains by 
covalent bonding along the field direction. 

Hydrogen atoms in strong magnetic field have been studied extensively We have 

also recently completed a study of the electronic structure of hydrogen molecules and chains 
in the strong field regime {b » 1) f^. However, in many of these studies, the center-of- 
mass (CM) effects of the proton motion have been neglected. In the case of the hydrogen 
atom, although significant effort has been devoted to calculating the energy levels of an 
electron in the static Coulomb potential of a fixed proton (infinite mass) to a high precision 
and for arbitrary magnetic field strength the two-body problem in strong magnetic 

field, including the effects of finite proton mass, has been studied in detail only recently 
(e.g., @Jl3l)- In this paper we focus on the hydrogen atom, but discuss those aspects which 
are important for application to molecules in very strong fields (we shall study molecular 



excitation levels in a later paper |TT[). 
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A free electron confined to the ground Landau level, the usual case for 6 >> 1, does 
not move perpendicular to the magnetic field. Such motion is necessarily accompanied by 
Landau excitations. When the electron (still in the Landau ground state) combines with a 
proton, the mobility of the neutral atom across the field depends on the ratio of the atomic 
excitation energy and the Landau excitation energy hujp = heB/{mpc) for the proton. As 
the typical Coulomb excitation is (In 6) in atomic units, it is convenient to define a critical 
field strength Bent via 

hcrit = — Inherit = 1.80 X 10^ B,rit = hritBo = 4.23 X 10^3 G. (1.2) 
rUe 

For B « Bo (terrestrial conditions) the Landau energies of both electron and proton are 
minor perturbations and one can construct wave packets that mimic the classical motion of 
a neutral atom across the field. In this case the internal structure of the atom or molecule is 
coupled to the center-of-mass motion only through a "motional Stark effect" (see Sec. II). In 
the opposite limit of i? >> B^rit » B^, the Landau energies are larger than the Coulomb 
excitation energy, and both electron and proton (in the atom) are in the Landau ground 
state at reasonable temperatures. In this case, quantum mechanics cannot mimic classical 
motion. The canonical field strengths of radio pulsars, with B slightly less than B^ru, present 
an intermediate case. However, at all field strengths one can introduce a pseudomomentum 
perpendicular to the field and, in principle, find the bound-state energy £^ as a function 
of Kx. One question of interest is the range of K± for which S increases linearly with K]_ 
(as it does for ordinary kinetic energy); but in any case S{K±) is needed to evaluate the 
Saha equation for the equilibrium between neutral and ionized hydrogen (Sec. V). Previous 
treatments [p!^ -[T5[| of the ionization equilibrium in strong magnetic fields did not properly 



take account of the non-trivial effects of atomic motion. States with large K_i_, where velocity 



decreases with increasing K± (see Ref. and our Eq. [3.29]) are of particular interest, 
especially for B larger than Bcrit- 

The separation of the center-of-mass (CM) motion of a bound state in the presence of a 
magnetic field has been studied previously based on the conserved pseudomomentum (e.g.. 



p!6|-pll). We briefly review and clarify this pseudomomentum scheme in Sec. 11. When 



B « Bcrit, perturbation calculations for hydrogen atom motion (e.g., |l9| , p2| , p3[ ) are valid 
over a wide range of K± values and lead to interesting phenomena such as bent trajectories 
p3| . Model atmospheres can be affected by details of the transverse motion (e.g., [p^ , p5[1 ). 



Some accurate numerical calculations for general K_i_ and for several values of B (but all 
below Bcrit) are now available but we concentrate on the B > Bcru regime in Sec. III- 



IV. A different approach to the two-body problem (for positronium atom) in the superstrong 
fleld regime 6 >> 1 has been developed in ||26|| . 

Our purpose in this paper is not to obtain accurate energy spectra of a moving hydrogen 
atom in certain limited regimes, as have been done in most of the papers mentioned above; 
rather, we seek complete (though approximate) solutions of the two-body problem for a 
wide range of parameter space, including B just below Bcru (a common case for neutron 
stars), but especially for B > Bcrit (in case B > 10^^ G exists in some neutron stars, as 
has been suggested recently p7|-p9|]). Our emphasis is on finding physically meaningful 
approximate fitting formulae for the atomic energy of the moving atom over all relevant 
values of Kj_ (Sec. Ill), in order to determine the equilibrium between neutral and ionized 
hydrogen (Sec. V). However, in practice we shall be interested mainly in the regimes where 
the thermal energy fc^T is much less than the ground-state binding energy of the atom, 
while the gas density is much smaller than the internal density of the atom, so that the 
neutral and ionized fractions are of the same order of magnitude. In such cases, we are most 
interested in two kinds of excited bound states: (i) Those with excitation energies up to a 
few atomic units (comparable to ksT but a small fraction of the binding energy) and (ii) 
states that are only barely bound (e.g., those with extremely large K_i_), for which one has to 
check whether phase-space factors make them unimportant relative to ionized hydrogen. For 
applications to molecules @Jll[] and multi-electron atoms |^ with B » Bcrit-, a controversy 
arises regarding the "coupling" of the electron's orbital quantum number with the Landau 
level of the proton (or nuclei). This is discussed in Sec. IV. 

Our calculations in this paper are based on non-relativistic quantum mechanics. For 



B > Brei = (hc/e^YBo = 4.414 X 10^^ G (note that Brei is close to Bcru only by coincidence), 
i.e., hue ^ frieC^, the transverse motion of the electron becomes relativistic. However, the 
relativistic correction to the atomic binding energy is small as long as the electron remains 



non- relativistic along field direction [^. Except as otherwise noted, we shall use atomic 
units (a.u.) throughout the paper, in which mass and length are expressed in units of the 
electron mass me and the Bohr radius = h"^/ (meC^) = 0.529 x 10^* cm, energy in units of 
2 Rydberg = e^/oo = 2 x 13.6 eV; field strength is in units of Bo (Eq. [1.1]) and temperature 
in units of 3.15 x 10^ K. 



II. SEPARATION OF CENTER-OF-MASS MOTION: THE 
PSEUDOMOMENTUM APPROACH 

To set the scene, we briefiy review the pseudomomentum approach to the two-body 
problem of a hydrogen atom in a strong magnetic field [|r^-|TU[. However, the physical 
meaning of the pseudomomentum of the atom needs some clarification. 



A. Pseudomomentum 

For a free particle of charge Cj and mass rrii in a constant magnetic field (assumed to 
be aligned in the z-axis), there are three momentum-like vectors: the canonical momentum 
P = — iV; the mechanical momentum 11 = P — CjA = mjV, where A is the vector potential 
and V is the velocity; and the pseudomentum (or the generalized momentum), as defined by 

K = n + eiBxr. (2.1) 

That K is a constant of motion can be easily seen from the classical equation of motion 
for the particle dll/dt = ei{dr/dt) x B. The parallel component Kz is simply the linear 
momentum, while the constancy of the perpendicular component is the result of the 
fact that the guiding center of the circular orbit of the particle does not change with time. 
The position vector Rc of this guiding center is related to by 
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Mathematically, the conservation of K is the result of the invariance of the Hamiltonian 
under a spatial translation plus a gauge transformation 0]. 

The existence of the integration constant or Rc implies infinite degeneracy of a 
given Landau energy level. We can use K to classify the eigenstates. However, since two 
components of do not commute, [K^^, Ky] = —iCiB, only one function of K;^., Ky can be 
diagonalized for stationary states. This means that the guiding center of the particle can 
not be specified accurately. If we use to classify the states, then the wavefunction has 
the well-known form e^'^^(j){y) where the function (f){y) is centered at yc = —Kx/{eiB), 
as can be inferred from Eq. (2.2). The Landau degeneracy in an area Ag = L"^ is thus given 
by {Lg/2TT) jdKx = {Lg/2iT)\Kx^g\ = Ag{\ei\B/2n), where we have used K^^g = —CiBLg. 
On the other hand, if we choose to diagonalize Kj_ = + Ky, we obtain the Landau 



wavefunction H^nm(r±) in cylindrical coordinates |^2|, where m is the "orbital" quantum 



number (denoted by s in some references). For the ground Landau level this is 

lVo™(r,) . = ^^J^ [Qj -M-^me). (2.3) 

The distance of the guiding center of the particle from the origin is given by 

= (2m +1)1/2^, m = 0,l,2,---, (2.4) 
where p is the cyclotron radius 

p = f^y^' = a„ f = b-'/^ (a.u.) = 2.57 x IQ-'^B;,'^' (cm). (2.5) 



eB J "\B 

The corresponding value of K± is given by = \ei\B{2m + 1). Note that Kj^ assumes 
discrete values since m is required to be an integer in order for the wavefunction to be single- 
valued. The degeneracy rrig of the Landau level in an area Ag = irRg is then determined by 
pmg — {2mgY/'^p = Rg, which again yields nig = Ag\ei\B /{2n). 

We also note that is related to the z-angular momentum Jz, as is evident from the 
g-jm6» fa^Q^Qj. in cylindrical wavefunction (Eq. [2.3]). In general, we can show that 



Jz = xpy - yp^ = ——{Kl - n^) = (m - n) — , (2.6) 
zeiiD ti 

where we have used n5_ = \ei\B{2n + 1), and n is the quantum number for the Landau 

excitations. 



B. Hydrogen Atom as a Two-Body Problem 

We now consider the electron-proton system. It is easy to show that even with the 
Coulomb interaction between the particles, the total pseudomomentum is a constant of 
motion 

K = Ki + K2, (2.7) 

where the subscripts 1,2 refer to electron (charge —1) and proton. Moreover, unlike the 
single particle case, here all components of K commute. Thus it is natural to separate the 
CM motion from the internal degree of freedom, using the vector pseudomomentum K as 
an explicit constant of motion (although we shall discuss an alternative set of basic states 
in Sec. IV). From Eq. (2.2) we have 

Rii" = Rci — Rc2 = — • (2.8) 

Thus we see that K_l is directly related to the separation Hk between the guiding center of 
the Landau orbit of the electron and that of the proton. 

Consider the energy eigenstate with a fixed K. Introduce the center-of-mass coordinate 
R = {meri+mpr2) /{me + nip) and the relative coordinate r = ri — r2. Writing the two-body 
wavefunction as 



^(R, r) = exp 



i(K + -B X r) ■ R 

, 2 ' . 



(2.9) 



so that \&(R, r) has a well-defined value of K, while exp[z(B x r) ■ R/2] is a gauge factor, 
the Schrodinger equation reduces to Q 



^The spin terms of the electron and the proton are not explictly included. However, it should be 
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H<P{r) = {Ho + H')(l){r) = £ 0(r), (2.10) 

with 

K"^ I f I \^ 1 1 
i/o = — + — P + -Bxr B-(rxp)--, (2.11) 

iy' = ^ + — (KxB)-r, (2.12) 

where p = —id/dr, M = rUf, + rUp ~ mp/rrie (a.u.), /i = rriemp/M ~ mg. We will make 
no distinction between M and rUp, nor /i and mg in our following calculations. Eqs. (2.10)- 
(2.12) have already been derived in Refs. |16-19|- Clearly, the CM motion is coupled to 



the internal motion through the second term in H', which has the form of the potential in 
a mot ion- induced electric field (K/M) x B. This term represents the so called "motional 
Stark effect" (although such a description is not exactly accurate, since K_|_/M does not 
correspond to the CM velocity |2^). For small K±, this effect can be considered by treating 



H' as a perturbation (Sec. Ill), but the eigenstates of Ho can in principle be used as a set of 
basic states for developing eigenstates of the full Hamiltonian for any value of Kj_. However, 
the following transformed version of the Hamiltonian is more convenient for large K±. 

Motivated by the fact that measures the separation of the guiding centers of the elec- 
tron and the proton, we can remove the "Stark term" by introducing a displaced coordinate 
r' = r — Hk, where Hk is given by Eq. (2.8). After a gauge transformation, with 

0(r) ^ exp ( t'^-^K^ . r) 0(r'), (2.13) 



2M 

the Hamiltonian becomes 

K"^ 1 / 1 \^ 1 1 

ff = — + — p' + -Bxr' B-(r'xp')-i r, (2-14) 

2M 2fiV 2 J nip ^ ^' r' + R^' ^ ' 



understood that for the ground Landau state, the zero-point Landau energy is exactly cancelled 
by the spin energy. Also, the abnormal intrinsic magnetic moment of the proton is neglected, since 
it does not play a role in our analysis. 
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where p' = —id/dT'. This expression has been obtained in Refs. [|T8],|T9|. We shall see in 
Sec. Ill that this alternative form of the Hamiltonian is useful in the regime where h is much 
larger than hcrUi defined in Eq. (1.2). 



III. APPROXIMATE SOLUTIONS AND FITTING FORMULAE 

A. Zeroth Order Solutions 

We consider first the Hamiltonian formulation in terms of Eqs. (2.11)-(2.12). For the 
zeroth order Hamiltonian the quantum numbers for the basic states are K^-, the number 
of nodes in the z-wavef unction z/, the electron Landau level integer n and the "orbital" 
quantum number m = 0,l,2,---. In this paper we only consider 6 >> 1 and thus restrict 
ourselves to n = 0. The energy eigenvalues of Ho for the = states can be written as 

£(o)(Kj = ^ + Era + m— , (3.1) 

where is the energy of a bound electron in the fixed Coulomb potential of an infinitely 
massive positive charge. The last term in Eq. (3.1) for m > 1 represents Landau energy 
excitations for the proton, but m is merely the "orbital" quantum number for the electron 
wavefunction and measures the relative z— angular momentum Jz = —m. Thus, there is no 
separate quantum number for the proton in this formulation. The "coupling" between the 
electron quantum number m and the proton Landau excitation mb/M in Eq. (3.1) results 
from the conservation of total pseudomomentum. The term Em has the form (e.g., ll|,|]) 

Em^-0.l6All, l^ = \n—^, (3.2) 

2m + 1 

where A is a coefficient which varies slowly with b and m (E.g., A ~ 1.01 — 1.3 for m = — 5 
when Bi2 = 1, and A ~ 1.02 — 1.04 for m = — 5 when B12 = 10). In most formulae below, 
we replace A by unity; the numerical values can be found in PJT^ . The atom has a cigar-like 
shape, with size ~ pm (cf. Eq. [2.4]) perpendicular to the field and ~ along the field 
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direction. Note that the correction due to the reduced mass /i could be easily incorporated 
by a simple scaling: 



E ~ 



-0.16A 



In 



(3.3) 



2m + 1 ^ 

However, this is a small correction (of order me/mp), and will be neglected hereafter. 

Equations (3.2)-(3.3) refer to the "tight-bound" states for which the number of nodes 
u of the z-wavefunction f{z) of the electron is zero. For z/ > 0, the energy eigenvalues are 
approximately given by |33| : 



E„ 



2vl 



1 - 



4pm 



z/i = 1,2,3, 



(3.4) 



for the odd states [y = 2ui — 1), and 
E - 1 



2ut 



= 1,2,3, 



(3.5) 



z/i \n{ao/pm 

for the even states {u = 2vi). The sizes of the wavefunctions are pm perpendicular to the 
field and ~ z/^ (a.u.) along the field. These states have much lower binding energies 
compared to the "tight-bound" states. 

We now consider the energies and eigenstates of the atom for finite Kj^. The two different 
Hamiltonian forms are discussed in Sec. III.B-C, and the general approximate expressions 
for the energies are then given in Sec. III.D. We focus on the "tight-bound" states only, since 
finite i^^ will make the weakly-bound z/ > states even less bound (although in Sec. V we 
will include an estimate of the statistical weight of these states in the partition function of 
bound atom). 



B. The Perturbation Hamiltonian Formalism 

For sufficiently small Kj^, we can use standard perturbation theory to calculate the 
correction of energy £ due to H' given by Eq. (2.12) (see also [^,^,^). Let K_l be along 
the ?/-axis, then the r-dependent part of H' is K^hx/M. We consider only n = and 
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u — 0, so the exact eigenstates of Hg + H' are superpositions of the ifo-eigenstates with 
m — 0,1,2, ■ ■ ■. The only non-zero matrix elements of x are of the form 

{Wm\x\W^+,) = j P, (3.6) 

and the energy differences of adjacent i?o-eigenstates are approximately given by 

A£r-^iJi.-4°'-0.32nn(H^) + A. (3.7) 

where I — Iq = \nb and the factor 0.32 is an approximation to a slowly varying function 
of m and b. We first consider the ground state m = 0. Using H' = K±bx/M + K']_/{2M) 
and Eq. (3.6) we note that perturbation theory is justified if is much smaller than a 
"perturbation limit" K^p defined as 

where ^ ~ 2.8 is a slowly varying function of b (e.g., ^ ~ 2 — 3 for B12 = 0.1 — 10^). For 
Kj^ « K±p, the energy Sq^^ to be added to Sq'^ in Eq. (3.1) is given by second order 
perturbation theory (plus a diagonal term) as 

where M± is the effective mass for the "transverse motion" of the atom, 

Thus the effective mass M± increases with increasing b. 

Similar calculations for the m > states yield — K]_/{2M±m), with 

m + 1 m 



M ^ b 



b/M + 0.16 - 0.16 b/M + 0.16 - 0.16 /, 



(3.11) 



A convenient (but approximate) expression for the effective mass is given by 



M^^ ~ M + Ui'im + 1)-^, (3.12) 
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where is of the same order of magnitude as ^, but different (by a factor of a few) for 
different m-states. The important feature in Eqs. (3.11)-(3.12) is that the effective mass is 
larger for the higher-m state. 

The quadratic form of the effective "transverse kinetic energy" in Eq. (3.9) is vahd only 
when it is much less than the "perturbation limit", reached when K± = K^p. Using Eq. (3.8) 
and the approximation in Eq. (3.10), this kinetic energy limit becomes 



For 6 << bcrit (even if 6 >> 1) this limit is large compared with 1 a.u., so that the quadratic 
perturbation energy (or transverse kinetic energy) in Eq. (3.9) is valid for the most important 
(low energy) states. Moreover, from Eq. (3.10) the effective mass is close to the actually 
proton mass M. For superstrong fields, b » bcru, on the other hand, the effective mass (for 
K± << K±p) is much larger than M and the perturbation formalism already breaks down 
when the transverse kinetic energy is only of order 1 a.u. 

At least for B » Bcrit-, we have to consider values of K± large enough so that the per- 
turbation treatment for the formalism in Eqs. (2.10)-(2.12) is unsuitable. For any magnetic 
field strength there are still eigenstates for arbitrarily large values of i^x, but the transformed 
Hamiltonian in Eqs. (2.13)-(2.14) is now more suitable for calculating the energy. 



As mentioned, Eq. (2.14) gives an alternative formulation for the Hamiltonian where 
K'j_/{2M) does not appear explicitly, but the displacement of the electron-proton guiding 
centers does with Rk = K±/b (in atomic units) ||T9| , p!0[| . We again focus only on the Landau 
ground state for the electron {n = 0), but in principle we must include all m- values (with 
proton Landau excitation energy mb/M) and mixing between these states. We consider 
first the approximation where we omit mixing, i.e., we use the diagonal matrix element of 
|r' + Rx|^^ for a fixed m- value, and restrict ourselves to u = 0, = and (2m + 1) << b, 




(3.13) 



C. The "Displaced Center" Formalism 
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so that = (2m + l)/b << 1. We can estimate the size of the atom along the 2;-axis 
and the energy for two different regimes of the values of Rk- 

(i) For Rk ^ Lz ^ I (but not necessarily Rk < pm) and with as a variational 
parameter we have 



1 1 J mb 



^m----\r.-,r^ + —^ (3.14) 



where the first term is the kinetic energy along the 2;-axis, and the second term is the potential 
energy of the electron. The logarithmic factor in Eq. (3.14) comes from an integration 
over the "cigar-shaped" electron cloud in the displaced Coulomb potential — l/|r' + Ric-]- 
Minimizing with respect to L^, we obtain 

The mixing between different m-states is unimportant when b » bcru- This can be seen 
from the order of magnitude estimate of the off-diagonal matrix element between m = and 
m = 1 states: 







1 



|r' + R 



K\ 



1 > ~ ( 



Rx r' 



(3.16) 



as compared with /S.£rn ~ {l + b/M) ~ l{l + b/bcrit)- For b ^ bent the mixing is non- negligible, 



especially when Rk ~ p << 1; some results are given in Ref. |jTO[- When 6 >> bcru-, the 
mixing can be neglected for all Rk- 

(ii) For Rk ^ 1, the Coulomb logarithm in Eq. (3.14) disappears, and we have (for 
m = 0): 

^» ~ If - (LfrW^- 

In the limit of Rk >> 1, minimization of £q with respect to yields 

Lz ~ rT « Rk. £o ~ = (3.18) 

Kk a_l 

Thus for K± ^ b, the atom is very weakly bound (|£^o| ^ !)• The limiting scaling relations 
in Eq. (3.18) have been identified in 
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We can calculate the energy eigenvalue more accurately. Substituting 0(r') = 
W^m(r^)/m(-2) into H(j){r') = S(j){r') with H given by Eq. (2.14), and averaging over the 
transverse direction, we obtain a one-dimensional Schrodinger equation 

1 

- ^^•^™(^) + ^K)fm{z) = £mfm{z), (3.19) 

where the averaged potential is given by 
Vmiz^Rx) = 





1 







WUr'^)j. (3.20) 
The function Vm{z,RK) can be evaluated using an integral representation (e.g., |1TT|| ) 

Vm{z, Rk) = -T dqexp l-— - g-^ 1 Jo{qRK/p)Lm i—j , (3.21) 
where Jo is the Bessel function of zeroth order and Lm is the Laguerre polynormial of order 



m 



0. We solve for E by integrating Eq. (3.19) numerically from z = oo to z = subject to 
the boundary conditions df /dz = at z = and / ^ as z ^ oo. The energy eigenvalue 
Sm as a function of Rk are shown in Figure 1 for B12 = 10 and 100. Similar numerical 
results have also been obtained in Ref. [|T^. For Rk ^ 1, we can fit the energy to a form 
similar to Eq. (3.15). The total energy of the m = state is then given by 

S,(,UK,,^§-0.w(ln^^)\ (3.22) 

and the atomic size ~ 1/| \ia{p^ + CRj^)\. Equation (3.22) reduces to Eq. (3.2) for K± = 0. 
From the numerical results we find C ^ 0.8. For Rk ^ 1, the binding energy of the atom is 
much smaller, and Eq. (3.22) should be replaced by the order of magnitude relation (3.18), 
while the actual numerical values of the energies are not important in practice (see Sec. V). 

In the B » Bcru regime, the maximum value K±p for the perturbation treatment of 
Sec. III.B to be vahd becomes K^p ~ b^/^ (cf. Eq. [3.8]). Thus for « 6^/^ Eq. (3.22) 
should be consistent with Eq. (3.9). Indeed, when Rk << p, Eq. (3.22) reduces to 

SoiK,, K^)^^- 0.16 /2 + 0.32 C^-Kl. (3.23) 

15 



The dependence on is again quadratic, and the corresponding effective transverse mass 
is ~ h/{2ACl) ~ 26// ~ ih/l, in agreement with Eq. (3.10). 

The m > states can be similarly calculated using Eqs. (3.19)-(3.21). Some numerical 
results are again shown in Figure 1. The energy can be expressed approximately as 

A-x) - § + ™^ - 0.16 (in -r^^ ' ■ (3.24) 

From the numerical results we again obtain Cm ~ 1- 

Comparison with the numerical results of Potekhin [0], who included the mixing of 
different m— states, indicates that Eqs. (3.22) and (3.24) are accurate to within ~ 30% 
in the relevant regime of i^_L (< h) when h > bcrit- The agreement becomes better as b 
increases. For smaller b, however, the perturbative results of Sec. 3.B should be adequate 
(see Sec. 3.D). 

Finally, if we consider sufficiently strong magnetic field so that not only b/M > I (or 
b ^ bcrit) but also b/M ^ l'^ is satisfied, then Eq. (3.24) implies that all the m > states are 
unbound, as has already been noted in Ref. [jl9|. However, this does not mean that there 
is no other bound state except a single non-degenerate m = state. Indeed, Eq. (3.22) 
indicates that there are many states for which the guiding centers of proton and electron 
are separated by a small Rk, and these states have similar energies compared to the ground 
state (m = 0, Kj_ = 0). In the pseudomomentum scheme discussed here, these states occupy 
a continuum i^_L-space. As we shall see in Sec. IV, these closely-packed energy levels can be 
made discrete if we use a different set of eigenstates. 



D. General Fitting Formulae 

Consider for the moment the cases with Kj_ << b, i.e., Rk << 1, but with no other 
restrictions on K_i_ or b. In Sec. III.B-C, we have obtained reasonably accurate ground-state 
(m = 0) energy of a hydrogen atom in two limiting regimes: (i) for B << Bcru, where 
Eqs. (3.9)-(3.10) are applicable up to adequately large values of K±, and (ii) for B » Bcru, 
where the energy is given by Eq. (3.22). We write the total ground state energy in the form 
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8o{K,, ir^) ~ ^ - 0.16 /2 + E^{K^), (3.25) 

and want to find a general fitting formula for the "transverse kinetic energy" E±{K±) with 
K± « b, but K_i_/b^^'^ = Rk/p otherwise arbitrary. With the inequality Rk « 1, the 
second term in Eq. (3.22) can be approximated by —0.16/^ + 0.32/ ln(l + CR]^/p^). We 
propose the following fitting formula 

^±(^±)^^ln(l + ^), r^0.64eiri^, (3.26) 

where is given by Eq. (3.10). Recall that the parameter ^ ^ 2 — 3 for B12 = 0.1 — 10^, 
and a typical number to use is ^ = 2.8. This formula reduces to Eq. (3.9) when K± « Kj_p 
and to the above approximation of Eq. (3.22) when B » Bcru- We expect Eq. (3.26) to 
be accurate to within 30% for K± « b. In the regime b >> bcrit and 6 << Kj_ « b^ 
(i.e., p « Rk « 1), Eq. (3.26) reduces to E± ^ 3.1 \n{Kl/b) (a.u.) >> 3 (a.u.), which is 
large compared with the thermal energy (temperature) T for the astrophysical applications 
of interest. For evaluating the integral over K± extending from zero to infinity in the atomic 
partition function (Sec. V), we shall advocate using Eq. (3.26) for all K± even though this 
expression tends to infinity, whereas the correct expression should approach the finite 
limit 0.16 for K± » b. The difference is appreciable only where E± » T (so that 
the Boltzman factor e~^^/'^ is very small) and our prescription amounts to "cutting off" 
the integral, i.e., omitting the states with Rk >> 1 from the integral. This omission is 
advantageous, since these states should be treated together with the ionized states, and 
both turn out to be unimportant (Sec. V). 



slightly more general fitting formula which closely resembles Eq. (3.22) is given by 

{l + Kl/Tf\ ' 

with r = 1 + Ml/{^b) = r/(0.64/M^). This expression can be applied to K±^ ^ b and gives the 
correct limiting result for K± >> b, but it is not convenient to use in practice. For the applications 
discussed in Sec. V, Eqs. (3.25)-(3.26) are adequate. 
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0.16 



Our fitting formulae for the energies of the m > states are less accurate. In the small- 
limit, the i^j_-dependent term in Eq. (3.24) reduces to the quadratic form K\/ (2Mx^), 

with the effective mass given by Eq. (3.12). Similar to Eqs. (3.25)-(3.26) we fit £rn to the 

analytical form: 

£^{K,, K^)c^^ + m^- 0.164 + E^^K^), (3.27) 



where 

E±m(Kj_) ~ In ( 1 + ^ ) , ^ 0.64e^(2m + 1)6 



U(2m + l)6j ' ^^-^^^ 
so that E±Q = E±, tq = T. Although can vary by a factor of a few for different values 
of B and m, equation (3.28) has the correct approximate functional dependence for a wide 
range of B and K±. 

As noted in the introduction, one cannot construct a wave packet to mimic the classical 
behavior of a moving atom when both electron and proton are confined to the ground Landau 
level. Therefore when B » Bcru » -Bo, the notion of "motion across the magnetic 
field" does not have a unique meaning. Nevertheless, one can ask two questions about the 
transverse pseudomomentum: (i) Is there an upper limit to K± beyond which there is no 
bound state? (ii) When a bound state exists, what is the value of Vj_ = dE/d'K±, the analog 
to the classical center-of-mass transverse velocity of the atom? From our discussion in this 
section, we have seen that there exist bound states for all values of B and K±, although the 
states with large K± are very weakly bound. Prom Eq. (3.26), we have 

(Note that for K± >> b, Eq. [3.18] should be used and wc have V± ~ i>/ K\ instead). Clearly, 
in general Vj_ is smaller than its field-free counterpart: V±_ ~ K^/M±_ (so that the effective 
mass description is valid) only for K^^ < K±p, and M± ~ M (classical behavior) only for 
b << bcrit- As K± oo, the center-of-mass velocity approaches zero. For a given magnetic 
field strength, the maximum V± is given by 

V±max ~ ~ ~ ^ b-^^^^ (3-30) 

2Mi Ml ' ^ ^ 



which occurs when ~ r^/^ ~ K±^p. For h ^ fccrii, the states with » K±p have not 
only small velocity but also large electron-proton separation » p. Nevertheless, these 
states are quite distinct from an electron-proton pair with separation Rk in fully ionized 
hydrogen, because the relative 2;-coordinate satisfies a bound-state wavefunction (localized, 
although with large scale- length and small binding energy). The partition function for these 
states is smaller than that for the ionized component, because the sum over the relative 
momentum in the z-direction, k^, is absent. 



IV. ALTERNATIVE APPROACH TO THE TWO-BODY PROBLEM 

The basic states used in Sec. Ill for the e-p two-body problem are explicit eigenstates 
of the transverse pseudomomentum K^. For K± = 0, this formulation has the advantages 
that the electron's "orbital" number m is a good quantum number, the wavefunction can 
be related to that for a fixed positive charge and the energy is given explicitly by Eq. (3.1). 
The K± = states with different values of m are orthogonal to each other and could be 



used to satisfy the Pauli principle for the electrons in a hydrogen molecule ||§|Tl[| or a multi- 
electron atom The last term mb/M = m{b/bcrit) In fecri* in Eq- (3.1) is unimportant when 
b << bcrit, but it is very large when b >> bcrit (and m > 1). The simplest wavefunction for 
the H2-ground state in this formulation would use one m = and one m = 1 electron |p, and 
the energy would include the positive term {b / bcrit) ^"^^ bcrit- One therefore might conclude 
that H2 has a positive energy relative to two H atoms when b » bcrit- In a future paper 



11| we will show that this is not the case. In a molecule, only the total pseudomomentum 
is conserved, not that for individual atoms. Thus for b >> bcrit one would have to perform 
an integral over of individual atoms to form a molecular eigenstate (Similar situation 
exists for multi-electron atoms). To our knowledge this complicated task has not yet been 
carried out. We propose here an alternative set of basic states which are not eigenstates 
of the pseudomomentum. This formulation is more complex for single H atom, but can 
be generalized more easily to molecule. Moreover, this formulation makes the definition of 
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eigenstates clearer (unlike the pseudomomentum approach, where m is a good quantum only 
for K± = 0), and therefore offers an intuitive understanding of the degeneracy of states. 

Let ri, T2 be the position vectors of electron and proton in a H atom. We introduce 
coordinates Z — {rUeZi + mpZ2)/M and z — zi — Z2i but stick to ri^ and r2_L. In this mixed 
coordinate system, the Hamiltonian for the electron-proton pair becomes H — Hg + V, with 
(cf. footnote 1): 

K"^ 1 
Ho = — + — +y (Pi± - CiAif, (4.1) 

" i=e,p ' 

where Aj = A(rj) = B x rj_|_/2, and the interaction potential is 

^""[^2 + (ri^_r,^)2]i/2- (4-2) 

We set Kz = Q without loss of generalality. 

Consider the transverse part of the wavef unction for a bound state. The Landau wave- 
functions of the electron and proton form a complete set 



nimi(ri±)VF'4„2(r2±), ni,n2,mi,m2 = 0, 1,2,--- (4.3) 

(^1,^2 specify the Landau excitations, mi,m2 are the "orbital" quantum numbers as in 
Eqs. [2. 3]- [2. 4]). In general, an eigenstate of H can be constructed as 

*(^,ri±,r2±) = Wn,mA^l±)W*^^^{v2A.)fn,m,n^m^{z), (4.4) 

ni,mi,n2 ,m2 

where we have restricted ourselves to z-wavefunctions without a node {v — Q). Substituting 
Eq. (4.4) into the Schrodinger equation {Ho-\-V)^ — and using the orthogonal relations 
for the functions W to average over the transverse direction, we obtain a set of coupled 
differential equations from which the functions /(-z)'s and the eigenvalue £ can be calculated, 
at least in principle. 

This set of equations is greatly simplified as a result of the conservation of total 2;-angular 
momentum J^. Form Eq. (2.6), we have 

= -(mi - ni) (m2 - 712). (4.5) 
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Indeed, since the basis function Wnirm{ri±)W*^^^{r2±) oc e-^('"i-"i)ei+i(m2-n2)e2^ ^hjle 
interaction potential V depend only on (6*1 — ^2), we readily see that only the states with 
the same Jz are coupled. 

We shall use the formahsm of this section only when B » Bcru, in which case the 
Landau energy of the proton is large compared to the atomic Coulomb energy, so that 
both Hi and n2 become good quantum numbers. For astrophysical applications we are then 
interested in the ground Landau levels, rii — n^. — 0. The energy eigenstate with a fixed 
^-angular momentum — rrit — m^. — mi can be constructed as 

*(-)(^,ri^,r2x) = EW^-i(ri±)W^;.+mi(r2±)/i7n^)- (4-6) 

mi 

The equations for the functions /'s are then given by 

- - E <i;(^)/iT^(^) - mi = 0, 1, 2, . . . (4.7) 

where 

<mi(^) = (^mi(m)W^;,+^,(r2±)|T-^|W^m;(rix)l^:,Wi(r2±)). (4.8) 
1 |ri — r2| 

The set of equations (4.7)-(4.8) essentially forms a differential-integral equation system (the 
sum over mi can be considered as an integration). Some mathematical formulae for evalu- 
ating the function G^^*l^, {z) in Eq. (4.8) are given in Appendix A. 

Since the states with different rrit are orthogonal, we can use the variational principle 
for each value of mt separately to find a rigorous upper bound to E^"^*\ To this end, 
we choose as a simple trial wavefunction the first term in Eq. (4.6), i.e., we include only 
mi = 0, m2 = rrit. Equation (4.7) then reduces to a single differential equation which is 
straightforward to solve numerically. For m^ = — 4, we find that the upper bound can be 
written in the form: 

£(»)<_o.i6^(„_i^y, (4.9) 

with C ~ 0.9 — 1.1 (depending on the values of rrit), and A is a coefficient close to unity (as 

1 /2 

in Eq. [3.2]). This expression is equivalent to Eq. (3.22), with m/ p playing the role of Rk, 
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but here rrit is an integer. The form in Eq. (4.9) is expected since rrit = m2 — mi measures the 
difference between the distances to the origin of the guiding center of the electron and that 
of the proton. Although the rrit > 1 state is not an exact i^'^-eigenstate, in a qualitative 
sense the separation between electron and proton increases with increasing rrit. Thus we 
expect 

which then naturally leads to the form in Eq. (4.9). The decrease from h in Eq. (3.2) to 
6/2.3 in Eq. (4.9) arises from the fact that both the electron and proton wavefunctions have 
finite spread around the same axis (unlike the usual pseudomomentum approach, where the 
relative coordinate is used). 

The actual energy eigenvalue £'(™-«) can be obtained by solving iteratively the series of 
equations in (4.7) using the standard shooting algorithm We have carried out the 



calculations for the rrit = and rrit = 1 states. Our numerical results are given in Figure 2 
for three different values of field strength B12 = 100, 1000 and 5000. Typically, more than 
10 terms in the sum of Eq. (4.6) are needed in order to attain convergence of the energy 
to within <, 1%. We find that the ground-state energy eigenvalue S^'^^ agrees with the 
standard value (Eq. [3.2] with m = 0) in the limit of B » Bcrit- This is expected because 
the rrit = ground state is also a K± = eigenstate (Note that K_|_i = — K_|_2 implies 
mi = m2 and hence rrit = 0). Also, we see that as B increases, the upper bound given 
by Eq. (4.9) becomes asymptotically closer to the actual energy £^"^*\ For general rrit, we 
cannot give a rigorous lower bound to E^"^*\ but it is approximately given by Eq. (4.9) with 
(6/2.3) replaced by b, so that the fractional uncertainty decreases as B increases. Note that 
although the contribution of the of mi > 1 terms in Eq. (4.6) to the correct energy becomes 
smaller as B increases, the contribution to the correct wavefunction is always non- negligible. 

It is instructive to consider the degeneracy of an energy eigenstate. Without the Coulomb 
interaction between electron and proton, there is a double Landau degeneracy in mi and m2. 
When the Coulomb interaction is included, the degeneracy in rut = m2 — mi is removed, but 
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a single degeneracy remains in mi, i.e., for a given rrit, the eigenfunction corresponding to 
is not unique: the eigenstate in Eq. (4.6) with real /mi's is presumably the state where 
the proton is centered at the origin; but there must be other states with the same energy, 
centered at different positions. We can demonstrate this degeneracy explicitly as follows. A 
Landau wavefunction centered at the origin of the coordinate can be expanded in terms of 
wavefunctions centered at some point Tq = ro± as 

Wnm{r^) = E e^^^'-^^'-Jn+m.n+m' (^] e'""""^/' W^nm' (r^ - T^), (4.11) 
m' \ " / 



'-n—m^n—m' 



where Fq = Xo + iyo = roe^^°, Kq is given by Tq = — Ko x B/i?^ (see Eq. [2.2]), /„ 
is the polynomial as defined in ||3^, and e*^" '"-^/^ is the gauge factor. We consider only the 
ground Landau level n = 0, and write 



(4.12) 



where C™' is a (complex) function of m, m' and Tq. The energy eigenstate in Eq. (4.6) can 
then be written as 



^('"*)(ri,r2) = e*Ko-(ri^-r2^ 



where we have defined 



)/2 



E 



W^m'i(rix - r,)iy;,+^,(r2x - ro)/i7"^*)(^) 



(4.13) 



mi 



(4.14) 



However, the term inside [■ ■ ■] in Eq. (4.13) is exactly an energy eigenstate [ri -To, r2- 
To) with = m[ based on the coordinate system centered at Tq. Since the state represented 
by the left-hand-side of Eq. (4.13) has a definite energy £^"^*\ while states with different 
have different energies, we must have m[ = nit in Eq. (4.13). Thus 



v^(™')(ri,r2) = e^K„.(rix-r.x)/2^(™0(ri - r„ - r„). 



(4.15) 



where we have added the subscripts "i, j" to indicate that there are many states associated 
with a given m^, i.e., the states with the same energy ^('"^^ is not unique. Clearly, the 
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degeneracy (per unit area) for a given rrit is B/{27i), i.e., a single Landau degeneracy (see 
the discussion following Eq. [2.5]). 

The above discussions demonstrate that there is a discrete set of states with rrit = 
1,2,3, ■ ■ ■, all having similar energies as the ground state {rrit = 0), and do not have any 
positive contribution rritb/M in their energies. This has important consequences for the 
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binding of hydrogen molecules in the B » Bcrit regime. In a forthcoming paper |TI 
shall use one rrit = and = 1 atom to construct the wavefunction for the H2 ground 
state, which also does not involve any Landau excitation of the proton. 

V. lONIZATION-RECOMBINATION EQUILIBRIUM 

A. Overview 

We now consider the ionization-recombination equilibrium of hydrogen atoms, e+p-v^ H, 
given by the generalized Saha equation in the presence of a strong magnetic field. Previous 
treatments of this problem (e.g., |12-|T5|) either ignored the coupling between the center- 



of-mass motion and the internal atomic structure, or did not have available our generalized 
formula for the "transverse kinetic energy" as a function of the pseudomomentum K±. 

Let T be the gas temperature in atomic units (about 3.16 x 10^ K), so that the Boltzmann 
constant fcs is set equal to unity, and Ug be the number density (also in a.u.) of protons 
(either free or bound) in the gas. We write Vg = irRg = l/ug and Ag = t^R^, so that a 
'Wigner-Seitz cylinder" of radius and length Rg contains one proton on the average. Some 
of the partition function integrals can be simplified if the density and temperature satisfy 
three inequalities: (i) The density is low in the sense that Rg = (vm^)"^/'^ is much larger 
than the largest dimension (i.e., the 2;-dimension) ~ of the ground-state atom; (ii) 
The temperature is much smaller than the ground-state binding energy \Eq\ = \E{H)\ ~ 
0.16/^ >> 1; (iii) The Coulomb attraction between a proton and an electron at typical 
separation Rg is of order R^^; we assume T » R^^ so that the "imperfect gas corrections" 
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are small. 

The Saha equation involves the bound-atom partition function Z{H) compared with the 
product Z[e)Z{p) of the two free-particle partition functions. Each of these systems has six 
discrete or "nominally continuous" quantum numbers. Two of these for the free e-p system 
refer to the z-motion, which can be represented by K^, the center-of-mass z-momentum and 
kzi the relative 2;-momentum. For the bound system (H atom), the i^^-partition function is 
identical, but instead of we have the quantum number v. In both systems the electron 
Landau excitations have the same quantum number n and energy nh, so that the n-partition 
functions are the same (In practice, 6 >> T, so we only need to consider the ground Landau 
state of electron and the n-partition function is essentially unity). For the bound system, 
the remaining three quantum numbers are m and the two Cartesian components of the 
transverse pseudomomentum Kj^. For the free e-p system the three quantum numbers are 
^2, the proton Landau level integer, and the two transverse parameters |i^±i| and |i^±2|- 



B. The Bound-State Partition Function of the H Atom 

We first consider only the ground state of the H atom, with m = u = 0. Using Eq. (3.25), 
the canonical partition function in the volume Vg — RgAg of a "Wigner-Seitz cylinder" can 
be written as 

Zm^R,{^)"\.,(^)z., (5.1) 

where the factor Rg{MT /27r)^/'^ results from the free center-of-mass z-motion and \E{H) \ ~ 
|£^o| — 0.16/^ is the ground-state binding energy. The partition function Z±_ associated with 
the "transverse motion" of the atom has the form 



Ag 



r^J-max 

/ 2TrK±dK± exp 
Jo 



T 



(5.2) 



(27r)2 

where E±{K±) is the generahzed "transverse kinetic energy". The upper limit K±jnax of the 
integral in Eq. (5.2) is determined by the condition Rk ^ Rg so that the pressure-ionized 
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states are excluded in the bound state partition function. It thus has the density-dependent 
form K±rnax ~ bRg. As discussed in Sec. III.D, E± is well approximated by Eq. (3.26) for 
K± up to ~ b. However, for K± >> 6 (or Rk >> 1) the correct expression for E± is close 
to 0.16 (almost independent of K±), whereas the approximate E± in Eq. (3.26) increases 
with K±. At very low density {Rk >> 1), there are highly excited states with K± between 
b and K±max ~ bRg, whose contribution to Z± is proportional to Agb^rig"^^^ exjp{— 0.161^ /T). 
However, these states cover only a narrow range of binding energies (of order Rg^ « T) 
and can be neglected compared with the ionized components, in view of the inequality 
0.16/^ » T. We can therefore omit these states entirely, but also make only a small error 
if we merely replace E± by the approximation in Eq. (3.26), which is an overestimate for 
b ^ Kx^ <^ bRg. This approximation has the advantage that the extension of the integral in 
Eq. (5.2) to infinity is not only finite but also small. We therefore get a convenient expression 
for Z_i_ by extending the integration to i^±maa; oo- We have 



Zl /■oo 
ZTT Jo 



A,^, (6.3) 



with 



Mi = M^(l-^^V\ (5.4) 



Thus Z± is proportional to M^, which is larger than M± (or equal to it when r/ {2M±T) » 
1, so that the effective mass approximation is valid throughout the regime of interest), and 
Mj_ is larger than the actual mass M. While the transverse motion is "slowed down" in 
the sense that dE±{K±)/dK± is smaller than the zero- field result, the K± ^ states still 
exist and their statistical weight oc M'^ is actually increased over the zero-field result by the 
strong magnetic field. 

We now consider the internal partition functions associated with the v > ^ and m > 
excited states, i.e., we write the total bound-state partition function as 

m) = y,{^^f'-^^^.{^^y(Ewm. (5.5) 
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Start with the quantum number v. The internal partition function relative to the ground 
state is 



Here is given by Eqs. (3.4)-(3.5) but can be approximated by = — l/(2z/^), and 
the factor of 2 comes from the near-degeneracy of the even and odd states. The sum in 
Eq. (5.6) is thus similar to that for a field- free atom except that the usual weight- factor 
2i/^ is missing as the atom is one-dimensional. Since \Eq\ — \Ey\ » T iox the temperature 
regime of interest, the individual term in the sum is very small, but the Zj,{H) would diverge 
if ^max were allowed to become infinite. As for the high Rydberg states in the field-free H 
atom, the divergence is avoided by including only states which fit inside the Wigner-Seitz 
cylinder. The size of the f >1 state is Lz ~ (a.u.), so we should choose Vmax ^ with 
energy Ei,^^^ ~ —l/{2Rg), i.e., we omit states with such extended wavefunctions that they 
would be pressure-ionized. With this prescription the sum in Eq. (5.6) becomes finite but 
density-dependent, i.e., it increases as n~^/® with decreasing Ug. Making use of the inequality 
T » 2\Ei,^^J, we split the sum into two parts: one extends from u — 1 to i/y, where 

= T; the other from vt to Vmax (which exceeds ut because of the inequality). The 
first part contains only a few terms, each with |£'o| — \Ei,\ » T, so they can be neglected 
(compared with unity, coming from the v — term). The second part represents the highly 
excited states and would contribute approximately 2 exp(— |£'o|/T)(7rnj,)~^/^, which could 
be large for very small density Ug. However, these states should be considered separately 
and compared with the ionized components: These states have negative energies of order 
{—1/Rg), while the ionized components have positive energies of order T. Because of the 
inequality T » Rg^, we neglect these states entirely in the rest of the paper and set 



The contribution of the m > states to the bound-state partition function can be 
considered in the same manner as that of the m = state discussed before. The internal 




(5.6) 
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partition function associated with the m-states is given by 



oo T\/fl 

m=0 ' 



(5.7) 



where 

M'^^ = M^m.{l-^^^^] , (5.8) 



so that M_Lo = Mj_ (cf. Eq. [5.4]), and M±rn, Tm are given in Eqs. (3.12) and (3.28). The sum 
in Eq. (5.7) simphfies in two extreme regimes of the field strength: For 6 << bcrit one may 
need to include a number of terms in the sum, but then M'j_^ is close to M; For b » bcrit 
the effective mass M'j_ is much larger than M, but since 6 >> MT one need to include only 
the m = ground state, and hence Zm{H) = 1. 

Finally, it is instructive to consider the partition function based on the alternative scheme 
discussed in Sec. IV. Using Eq. (4.9) (and relacing 6/2.3 by b; see the discussion following 
Eq. [4.10]), we see that Z{H) can still be written as Eq. (5.1), but the transverse partition 
function Z± is now given by 

A„b r"^t,; 



^ ,max 

mt=0 



l\n{l + C'mt) 



271 Jo 



/ dmt{l + C'mt)~''^, (5.9) 

JO 



where the factor Agb/2Ti comes from the Landau degeneracy of the mj-th state, and 
^t,max ~ R^b. Clearly, Eq. (5.9) has the same form as Eq. (5.3) in the limit of b/M >> I, 
demonstrating the equivalence of the two energy level schemes of Sec. II-III and Sec. IV. 



^The factor (1 + e~^^^'^'^) in Eq. (5.7) results from the proton spin term, which is not explicitly 
included in our calculations (see footnote 1). Note that when the proton spin term (but not the 
abnormal magnetic moment) is taken into account in Eq. (2.11) or (2.13), the energy of the m-th 
atomic state is given by £m ^ K^/{2M) - 0.16 4 + E±rn{K±) + mb/M + (1 + a^)b/{2M), where 
the proton spin cr^ = ±1 (compare with Eq. [3.27]). 
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C. Saha Equation 



The partition function of the free electrons in the ground Landau level is given by 

T\"= / . b 



where the factor i?g(T/27r)^/^ represents the free z-motion and the factor Agb/{2Tc) is the 
Landau degeneracy. For the free protons (see footnote 1), we have 

\ 



■MT\^/^ b 



exp 



2n'j 2n 



tanh 



MT , 

b " ^ 



2MT 



(5.11) 



where the sum extends over all Landau levels of the proton, and is the spin degeneracy: 
= 1 for n2 = and = 2 for n2 > 0. Given Z{e), Z{p) and Z{H), the ionization- 
recombination equilibrium can be obtained using the condition fi{e) + fi{p) = ^{H) for 
the chemical potentials. In the density and temperature regimes of interest, with T « 
Kip/(2M^), we have 

X{H) fbY\^ / 2MxT\-V^\'^' .fb\ f\E{H)\\ , 

^%(^) ^^(1-^) (^) tanh(^)exp(^^)^„(if), 

(5.12) 



XpXe \27iJ V T J \27rJ \2MT ) \ T J 



where X{H) = n{H)/ng, Xp = Up/ug, X^. = rie/ng are the number density fraction of 
different species, M± = M + ^b/l (with ^ ^ 2.8), and Zm{H) is given by Eq. (5.7). This is 
the generalized Saha equation in the presence of a superstrong magnetic field. More details 
on the applications of this result to neutron star atmospheres will be presented elsewhere 



VI. SUMMARY 



The effects of center-of-mass motion of neutral hydrogen atom in a strong magnetic 
field are rather intricate, mainly due to the high degree of degeneracy associated with the 
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quantum states. Using the usual pseudomomentum scheme (Sec. II-III), we have obtained 
approximate solutions for the energy of the atom as a function of the field strength and 
conserved pseudomomentum for a wide range of parameter regimes. In particular, we have 
considered field strengths B ^ Bcru ~ 4 x 10^'^ G, when the Landau excitation energy of 
proton is considerable. States with large transverse pseudomomentum have small binding 
energy and transverse velocity, but are nevertheless quite distinct from fully ionized states. 
We have concentrated on convenient analytic fitting formulae which give at least a reasonable 
approximation over various parameter regimes (see particularly Eq. (3.26) for the "transverse 
kinetic energy"). Since there may be neutron star atmospheres with B ~ 10^^—10^^ G, we are 
particularly interested in the cases with B » Bcru, where the proton Landau energy hup is 
very large. By considering an alternative scheme to the usual pseudomomentum formulation 
(Sec. IV), we have shown that there are atomic states with orbital wavefunctions orthogonal 
to that of the ground state, but without any Landau excitation appearing in their energies. 

We have also derived the generalized Saha equation for the equilibrium between neutral 
hydrogen atoms and the ionized component. We focused on the cases of astrophysical 
interest, where the density is relatively low and the thermal energy ksT is small compared 
to the atom's ground state binding energy. Although the maximum transverse velocity 
of bound atoms is small in strong magnetic fields, the statistical weight due to transverse 
motion is actually increased by the strong fields, not decreased (Sec. V). The statistical 
weight of highly excited bound states is smaller than that of the fully ionized component. 
Our results are important for determining the physical conditions of magnetic neutron star 
atmospheres as well as the soft X-ray (or EUV) radiation spectra from them. Some of these 
issues will be studied in a future paper |^ . 
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APPENDIX A: 

In this Appendix, we derive some mathematical formulae needed for evaluating the func- 
tion G^'jj/ (z) defined by Eq. (4.8). (Here the lengthscale is in units of the cyclotron radius 
p.) Using the identity 

1 1 fd^q 



we write Eq. (4.8) as 



27r2 



{rrit + mi\e^^^'^^^\mt + m'^)(m'^|e" 



(Al) 



(A2) 



Using the general result |^ for the matrix element {m'\e^'^^'^^\m) , and integrating over dg^, 
we obtian (assuming mi > m'l without loss of generality): 



2 \ mi — m-^ 



{mt + m'^)\m'^\ ^m^^m'J q^\ ^m^-^'J q^\ 



where is the Laguerre polynomial of order n [3^|. We now define constant coefficients 



(7^™''^(mi, m'^) via 



2, 



{mt +m'i)\ m'i\ m^-m.'JX\ m^-m'JX 

\\ {mt + mi)\mi\ "^*+k1v2; ™'i V2 



mi+m'j^+mt 

5: g^r\^i,^'i)Ln{x). (A4) 

n=0 



Using the relation 



dx- 



V2n\Jo ^^(x + ;2V2)V2' 
where the second equality can be used to evaluate the function V"„(z), we obtain 



(A5) 



(A6) 



n=0 



We calculate the coefficients (7^™*) (mi, m[) using the orthogonal relations of Laguerre poly- 
nomials. We can identify two special cases: (i) When mt = 0, we have G^^^^,^{z) = Em-^m[{z) 
and gl^\mi,m[) = e„(mi,m'^), where Emyrn[iz) and e„(mi,m'J are given in [||; (ii) When 
mi = m[ = 0, we have gI^''\z) = Domti^) and (7^™"'^ (0,0) = dn{0,mt), where Djnim[{z) and 
dn{mi,m[) are again given in f^. 



31 



REFERENCES 

[1] M. Ruderman, in Physics of Dense Matter (I.A.U. Symp. No. 53), Edited by C. J. 
Hansen (Dordrecht: Holland, 1974). 

[2] H. Ruder et al. Atoms in Strong Magnetic Fields (Springer- Verlag, 1994). 

[3] Early study of essentially the same problem dealt with excitons in semiconductors, 
where the strong field condition can be mimicked in virtue of the small effective mass 
and the presence of a dielectric medium. See R. J. Elliot and R. Loudon, J. Phys. Chem. 
Sohds 15, 196 (1960); H. Hasegawa and R. E. Howard, ibid. 21, 179 (1961). 

[4] For calculations in the superstrong field regime {B » 10^ G), see R. Cohen, J. Loden- 
quai, and M. Ruderman, Phys. Rev. Lett. 25, 467 (1970); V. Canuto and D. C. Kelley, 
Astrophys. Space Sci. 17, 277 (1972); J. Simola and J. Virtamo, J. Phys. B: At. Mol. 
Phys. 11, 3309 (1978), and references therein. 

[5] For a review of the intermediate field regime {B ^ 10^ G), see R. H. Garstang, Rep. 
Prog. Phys. 40, 105 (1977). 

[6] Hydrogen atoms in arbitrary magnetic fields were first studied by W. Rosner, G. Wun- 
ner, H. Herold and H. Ruder, J. Phys. B: At. Mol. Phys. 17, 29 (1984). 

[7] For high-precision calculations of the atomic binding energy, including the relativistic 
effects, see S. P. Goldman and Z. Chen, Phys. Rev. Lett. 67, 1403 (1991), and references 
therein. 

[8] D. Lai, E. E. Salpeter and S. L. Shapiro, Phys. Rev. A 45, 4832 (1992). 

[9] M. Vincke, M. Le Dourneuf and D. Baye, J. Phys. B: At. Mol. Phys. 25, 2787 (1992). 
[10] A. Y. Potekhin, J. Phys. B: At. Mol. Phys. 27, 1073 (1994). 
[11] D. Lai and E. E. Salpeter, Phys. Rev. A, in press (1995). 
[12] V. K. Khersonskn, Sov. Astron. 31, 225 (1987). 

32 



[13] M. C. Miller, Mon. Not. Roy. Astron. Soc. 255, 129 (1992). 

[14] J. Ventura, H. Herold, H. Ruder and F. Geyer, Astron. & Astrophys., 261, 235 (1992). 

[15] G. G. Pavlov, Y. A. Shibanov, V. E. Zavlin and R. D. Meyer, in Lives of Neutron Stars, 
ed. A. Alpar, U. Kiziloglu and J. van Paradijs (Kluwer, Dordrecht, 1994). 

[16] W. E. Lamb, Jr., Phys. Rev. 85, 259 (1952). 

[17] L. P. Gor'kov and I. E. Dzyaloshinskii, Sov. Phys. JETP 26, 449 (1968). 

[18] J. E. Avron, I. B. Herbst and B. Simon, Ann. Phys. (NY) 114, 431 (1978). 

[19] H. Herold, H. Ruder and G. Wunner, J. Phys. B: At. Mol. Phys. 14, 751 (1981). 

[20] B. R. Johnson, J. O. Hirschfelder, and K.-H. Yang, Rev. Mod. Phys. 55, 109 (1983). 

[21] Born-Oppenheimer approximation in strong magnetic field is discussed in P. Schmelcher, 
L. S. Gederbaum and H.-D. Meyer, Phys. Rev. A., 38, 6066 (1988). 

[22] M. Vincke and D. Baye, J. Phys. B: At. Mol. Phys. 21, 2407 (1988); D. Baye and M. 
Vincke, ibid, 23, 2467 (1990). 

[23] G. G. Pavlov and P. Meszaros, Astrophys. J. 416, 752 (1993). 

[24] G. G. Pavlov et al. Astron. & Astrophys. 289, 837 (1994). 

[25] J. Ventura, H. Herold, & Kopidakis, in Lives of Neutron Stars, ed. A. Alpar, U. Kiziloglu 

and J. van Paradijs (Kluwer, Dordrecht, 1994). 

[26] A. E. Shabad and V. V. Usov, Astrophys. Space Sci. 128, 377 (1986). 

[27] B. Paczyhski, Acta Astron., 42, 145 (1992). 

[28] R. C. Duncan and C. Thompson, Astrophys. J., 392, L9 (1993). 

[29] D. Chakrabarty, et al., in Proc. of the Second Compton Symposium, Edited by C. E. 
Fichtel, N. Gehrels and J. P. Norris (AIP Press: New York, 1994). 



33 



[30] D. Neuhauser, S. E. Koonin, and K. Langanke, Phys. Rev. A 36, 4163 (1987). 
[31] C. Angelie and C. Deutch, Phys. Lett. 67A, 353 (1978). 

[32] L. D. Landau and E. M. Lifshitz, Quantum Mechanics (Oxford: Pergamon Press, 1977). 
[33] L. K. Haines and D. H. Roberts, Am. J. Phys. 37, 1145 (1969). 

[34] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (New York: Dover) 
(1972). 

[35] W. H. Press, B. P. Flannery, S. A. Teukolsky and W. T. Vetterhng, Numerical Recipes: 
The Art of Scientific Computing (Cambridge Univ. Press, 1987). 

[36] D. Lai and E. E. Salpeter, Astrophys. J., in preparation (1995). 

[37] J. Virtamo and P. Jauho, II Nuovo Cimento, 26B, 537 (1975); A. A. Sokolov and I. M. 
Ternov, Synchrotron Radiation (Oxford: Pergamon Press) (1968). 



34 



FIGURES 

FIG. 1. Numerical results for the energies of hydrogen atom in the m = state (solid lines) 
and m = 1 state (dotted lines) as a function of Rk, obtained by solving Eq. (3.19). The upper 
curves are for Bu = 10, and the lower curves for Bu = 100. 

FIG. 2. The energy eigenvalue f ("^*) calculated from Eq. (4.7) for the mt = (upper panel) 
and mt = 1 (lower panel) states. The ratio £^'^*^/Eo, where Ep = — 0.16^Z^ is the ground-state 
energy as given by Eq. (3.2), is plotted against mimax, the maximum values of mi in the sum of 
Eq. (4.6). The filled triangles are for B12 = 100, the open circles for B12 = 1000, and the filled 
circles for Bu = 5000. 
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